Solitons with Cubic and Quintic Nonlinearities Modulated in Space and Time 
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This work deals with soliton solutions of the nonlinear Schroedinger equation with cubic and quintic non- 
linearities. We extend the procedure put forward in a recent Letter [J. Belmonte-Beitia et al., Phys. Rev. Lett. 
100, 164102 (2008)], and we solve the equation in the presence of linear background, and cubic and quintic 
interactions which are modulated in space and time. As a result, we show how a simple parameter can be used 
to generate brightlike or darklike localized nonlinear waves which oscillate in several distinct ways, driven by 
the space and time dependence of the parameters that control the trapping potential, and the cubic and quintic 
nonlinearities. 
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The nonlinear Schroedinger equation (NLSE) has been 
studied in a diversity of situations. The nonlinear interactions 
are usually of cubic nature, but there are systems which en- 
gender cubic and quintic (CQ) nonlinearities. The case of CQ 
nonlinearities opens new possibilities, and there are many in- 
teresting applications, specially in nonlinear optics, in fibers 
where the CQ nonlinearities may be used for instance to de- 
scribe pulse propagation in double-doped optical fibers, when 
the type of dopant varies along the fiber, with the value and 
sign of the cubic and quintic parameters that control the non- 
linearities being adjusted by properly choosing the character- 
istics of the two dopants [1]. There are other applications, and 
here we quote the Bose-Einstein condensates (BEC), where 
a diversity of nonlinearities may appear driven by controlled 
optical interactions 10]. 

A particularly important scenario for solitons concerns the 
NLSE in the case of a single spatial dimension with a linear 
term and cubic nonlinearity. However, we can change the lin- 
ear parameter to make it to trap the system. In this case it 
is named the Gross-Pitaevskii equation | 3], and the linear pa- 
rameter V is now modified to v{x), with explicit dependence 
on the spatial coordinate x. Usually, the potential v{x) is a 
background potential, and has the form of the harmonic po- 
tential, with the purpose of trapping the system in a finite re- 
gion in space, or some spatially periodically oscillating pat- 
tern, if the purpose is to entrap the system into a periodic lat- 
tice. This equation has gained further importance recently, 
mainly because of its direct application to the study of BECs 
Jil], in fibers and in photonic crystals and other periodic sys- 
tems Other applications include the study of nonlinear 
tunneling of spatial and temporal optical solitons in optical 
organic materials of practical use to ultrafast photonic tech- 
nologies 1 7]. 

The pioneer work by Serkin and Hasegawa (SH) fs*] has 
introduced an interesting procedure to deal with nonau- 
tonomous NLSE. It is based on a similarity transformation, 
which transforms nonautonomous NLSE equation into sta- 
tionary equation which is easier to solve. The SH procedure 
has inspired several investigations ||9|]. In particular, in the 
recent work by Belmonte-Beitia, Perez-Garcia, Vekslerchik, 



and Konotop (BPVK) d one deals with the equation 



dx'^ 



(1) 



Here v{x^t) and g{x^t) may now vary, being functions with 
both space and time dependence, the first one being the trap- 
ping potential, and the second one describing the cubic nonlin- 
earity. In the present work we shall follow the BPVK proce- 
dure, which is well explained in [10] and so we do not review 
it here. Instead, we extend the method to another problem, 
focusing our attention on the CQ nonlinear equation 



dt 



dx'^ 



-^x,t)^^gix,t)\^H^gix,t)m''^, (2) 



where v{x^t) is the trapping potential, and gs{x,t) and 
g^{x^t) control the nonlinear cubic and quintic interactions, 
respectively. The validity of the one dimensional approach is 
given explicitly in [10], and the quintic term which we have 
added in Q may simulate three-body collisions and/or devi- 
ation of the trapped condensate from the one-dimensionality 
LI 1]. The trapping potential and nonlinearities to be used be- 
low are typical of BECs, and the results obtained may stimu- 
late new experiments in the field. 

In this work, we provide a way of making the BPVK proce- 
dure work in this scenario, involving the presence of CQ non- 
linearities. The result allows one to obtain explicit solutions 
for some specific choices of parameters, leading us to interest- 
ing localized solutions of the bright or dark type, depending 
on the vanishing or not of the eigenvalue of the associated sta- 
tionary nonlinear equation, as we explain below. 

The idea is to write the solution of ^ as 



*(x,t)=/,(x,t)e^^(->*)$(C(x,t)), 
in order to rewrite © into the stationary equation 



(3) 



-+G3i*(c)r*(c)+G5i$(c)r*(c), (4) 



where ji is the eigenvalue of the nonlinear equation above, and 
G3 and G5 are real constants which control the CQ nonlinear- 
ities. The substitution of © into leads to ©, but now we 
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have to have 

dp d{p^{d4>/dx)) 



dx 



dt dx dx 



= 0, (5a) 

d{p'^{dC/dx)) 



dx 



= 0. (5b) 



We can introduce a new function i{x^t) such that C{x^t) = 
F{^{x^t)). In this case, we write S,{x^t) = ^{t)x -\- S{t). This 
choice is interesting because it allows to determine the width 
of the localized solution in the form l/j{t), and its center of 
mass position as —S{t)/'y{t). With the Eqs. (0 we can obtain 
the new set of equations 



{d^/dt) 2 {d5/dt) 



47 



27 



(6a) 

x + a{t), (6b) 



, ^, id^p d<t> fd<i>Y 
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(6d) 



g,{x,t) = G,G^^p-^gs{x,t)=G,j^p-^ (6e) 

where a{t) is an arbitrary function of time. The choice of F{^) 
is to be done in a way such that we get finite energy solutions 
of the CQ nonlinear equation. 

This is the general procedure, and we note from (l6el) that 
if we take the limit G5 ^ we get to gsixA) — » and 
this leads us back to the case set forward in llldl . This result 
is robust. In fact, we have shown that the procedure can be 
extended to other nonlinearities: if in Q we add the new term 
^7(x,t)|^|^^, etc, then in © we should add Grl^l^^, etc, 
and after ([6eb we should include the new expression ^7(x, t)= 
G7G^Y^g,{x,t)=GrG^Y^gs{x,t)=Grj^p-'^ and so 
on. We will deal with the more general case in another work 
I12I1 . and below we focus attention on the important case of 
CQ nonlinearities. 

Up the here, we have been very general. However, to il- 
lustrate the procedure with examples of interest, let us focus 
our attention into the case of specific nonlinearities. This is 
a nontrivial task, but we know that in the case of BECs with 
controlled optical interactions, we can have a diversity of non- 
linearities mSlllillSlTlli&llSIIil]. Thus, in the present 
work we suppose that the cubic nonlinearity is given explicitly 
by 



(7) 



where b is real parameter which controls the behavior of the 
cubic nonlinearity - we will use 6 = 8 where required. With 
this choice we have 



The quintic nonlinearity is now controlled by 



(8) 



and the potential given by Eq. (l6ct has to have the form 

v{x, t) =uj^x^+hx+h - pi^Gf/^'e^^"'^"" , (10) 
and cj, /i, and /2, are time-dependent functions such that 



u,\t) = 7" + 
/^(t) =72 +72^2. 
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(11a) 
(lib) 
(11c) 



We define x = I/7 in order to rewrite (II lab in the form of 
the Ermakov-Pinney equation iflolflS l 



df^ 



This equation has solution for 7 in the form 

^{t)=[2yl{t)^2yl{t)/W^Y^^\ 



(12) 



(13) 



Here W is the Wronskian of the two linearly independent so- 
lutions yi e y2 of the Mathieu equation {dP' y / dt^)^^uj'^ {t)y = 
0. See Ref. [13] for more details on this. 
We now choose 



uj'^{t) = 1 -\- £ cos(a;o^), 



(14) 



in order to obtain analytical solutions. The above results are 
general and can be used to investigate explicit examples, as 
we consider below. Before searching for explicit solutions, 
however, let us remark that investigations concerning specific 
criteria for adiabaticity of nonlinear wave equations and of the 
soliton solutions to such equations, and issues related to para- 
metric amplification of elementary excitations due to the pe- 
riodic modulation of the trapping potential and nonlinearities 
were studied before in and Ill5n . respectively. 

An important issue which appears in the case of CQ non- 
linearities is that we can solve Eq. ^ for bright and dark soli- 
tons, depending of the values of the parameters that control 
the nonlinearities in a specific system. As the first example, 
let us consider the simpler case in which the stationary Eq. ^ 
has vanishing eigenvalue, p = 0. In this case, if we choose 
the cubic and quintic parameters as G3 = 2 and G5 = —3, we 
can write Eq. ^ in the form 



2<I>^(C) -3$^(C). 



(15) 



The solution is given by ^{() = l/y^l + C^. It has the 
bell shape form, and it is of the bright type, as we plot in 
Fig. 1(a). To have the correct boundary condition, we have 
to set C ±(X) as X ±00, and so we have to choose 
C{x^t) = F(^(x,t)) properly in order to make the solution 
to behave according to the required boundary conditions. Ac- 



(9) cording to Eqs. Q, we choose F(^) = G3 / d^e' 
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FIG. 1: (Color on line) Plots of $(C) for the bright soliton (a) (blue, 
solid line) and the sech(C) (red, dashed line), and for the dark soliton 
(b) (blue, solid line) and the tanh(C) (red, dashed line). 
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FIG. 3: (Color on line) Plots of |^(x,t)|^ for the resonant bright 
soliton for e = 0.5 and uj^ — \ (a), and of the breathing bright 
soliton, for e = (b). Initial data for Eq. ([T2l ) are x(0) = \fi and 
dx/dt{0) = 0. 



We note that the bright solution is localized, but it is thicker 
then the hyperbolic secant, the standard bright soliton, which 
is also plotted in Fig. 1(a), for comparison. 

As we know, the form found above for the bright soliton can 
also be obtained in models of relativistic scalar fields, but this 
will be reported in another work [16] . We will get into this fol- 
lowing the lines of [17], which suggests several distinct ways 
of modifying the bell shaped form of brightlike solutions, an 
issue of direct interest for practical applications, as we will 
illustrate in [^]. 

To get to the darklike solution, let us now choose a nonva- 
nishing eigenvalue, /i 7^ 0. In order to get explicitly simple 
analytical solution we consider, for instance, the case /i = 3, 
G3 = 6, and G5 = — 3. In this case, Eq. dH) changes to 



= -3^(0+6^^(0-3^^(0. 



(16) 



The solution is given by $(0 = C/ \/l + C^- It the form 
of a kink, as we show in Fig. 1(b), and so |^(C)P is now of 
the dark type. It is interesting to note that this form of solution 
has also appeared in high energy physics, as reported before 
in II8]. It is thicker then the standard darklike soliton which 
is described by the hyperbolic tangent, which we also plot in 
Fig. 1(b), for comparison. 

Equation ([T5]) leads to the brightlike solutions. In this case 
the wave function which solves Eq. © acquires the form 



^(x,t) 



2l/6^1/2gi0g-^V66' 



(17) 





FIG. 2: (Color on line) Plots of the trapping potential v(x,t) given by 
JlSal ) for the resonant brightlike case (a) and by ( I20al ) for the resonant 
darklike case (b), in the range — 3 < x < 3 for the time evolution in 
the interval < t < 20. 



where (j) = (j){x^ t) is real, obtained via the Eq. ([6bb . 

A specific solution can be constructed, which leads to res- 
onant solitons. It is given by (5(t) = 0, /i = and a{t) = 
J ^'^{t)dt. In this case we get 



v{x^ t) 



-3(2) 



(18a) 
(18b) 



with ^3(x, t) given by d?]). In Fig. 2(a) we plot the trapping 
potential (II Sab in the standard form, to show how it behaves 
as a function of space and time. This is to be compared with 
the trapping potential of the other case, for /i 7^ 0, where the 
dark solitons appear. 

Solutions of the brightlike form nicely appear when one ad- 
equately chooses both e and uoq in Eq. ([T4l) . In Fig. 3, we show 
the resonant and breathing solitons for the specific choices of 
parameters e = 0.5 and cjq = 1, and e = 0, respectively. We 
can also get quasiperiodic and moving solitons, but we will 
leave this to the longer work lll2n . 

As we have already seen, the Eq. ([T6l) gives the darklike 
solutions. In this case, the Eq. ^ takes the form 



If we take 5(t) = 0, /i 7^ e a{t) = J -f^{t)dt, we get 
g,{x^t) = -U-'/'e'^'/''\ 



(19) 



(20a) 
(20b) 



with (x, t) given by (|7]). In Fig. 2(b) we plot the trapping 
potential (I20al) in the standard form, for the resonant solution, 
to show how it behaves as a function of space and time. Here 
we see that for the darklike soliton, the trapping potential has 
a different behavior, if compared with the case of brightlike 
soliton - see Fig. 2. The presence of a nonvanishing fi mod- 
ifies both the form and deepness of the trapping potential in 
the darklike case. 

In the darklike case, we also have several soliton solutions. 
The resonant dark soliton appears with an appropriate choice 
for e and uq in Eq. ([T4l) , as we show in Fig. 4. We can also 
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we have shown how to find the solutions with the resonant, 
breathing, quasiperiodic and moving behavior. A more de- 



FIG. 4: (Color on line) Plots of t) |^ for the resonant dark soH- 
ton for e = 0.5 and wo = 1 (a) and of the breathing dark soliton, for 
e = (b). Initial data for ^ as in Fig. 3. 



have the breathing soliton. It is obtained taking e = in 
([T4l) . In this case we have uj{t) = 1, and this leads us with 
j{t) = 2/v'l + 15cos2(2t). The t)p for darklike soli- 
ton of the breathing type is also shown in Fig. 4. If we take 
coq = >/2/2 in (O, we can use the above ^{t) to get to 
quasiperiodic solitons, as we show explicitly in Fig. 5. 

Up to here we have chosen S{t) = 0, to fix the center of 
mass of the solution; however, if we choose S{t)y^O we can 
make it to move. For simplicity, we consider the case e = 
in (O, and take j{t) as above. Thus, we get to fi{t) = in 
(lllbb and /2(t) = in (lllcb with the introduction of r(t) = 
/ 7^(t)(it, which leads to = cos(2r(t)) and a(t) =r(t) + 
(1/4) sin(4r(t)). With this, we get to the case of a moving 
dark soliton, as we show explicitly in Fig. 6. 
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FIG. 5: (Color on line) Plot of |^(x, t)|^ for the quasiperiodic dark 
soliton with e = 0.5 and cuo = a/2/2 in f]A[ (a). The quasiperiodic 
behavior is better seen from the soliton profile at panel (b). Initial 
data for ^ as in Fig. 3. 

In this work we have extended the procedure set forward 
in the recent Letter [10] to the case of cubic and quintic non- 
linearities. As we have shown explicitly, the procedure is ro- 
bust and can be generalized to higher (odd) nonlinearities very 
naturally, algorithmically. This result opens a new route for 
bright and dark solitons, with a diversity of possibilities of 
practical applications. 

To strengthen the result, we have shown explicitly that 
the system can support bright or dark solitons, depending on 
the way we make the eigenvalue of the stationary nonlinear 
Eq. d?]) to vanish, in the case of bright soliton, or not, in the 
case of dark soliton. In the specific case of dark solitons. 




FIG. 6: (Color on line) Plot of |^(x,t)|2 for the moving breathing 
dark soliton (a). The center of mass motion is better seen from the 
profile at panel (b). Initial data for ([T2l ) as in Fig. 3. 



tailed discussion is under preparation ||12|] . 
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